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ABSTRACT 

Context. Photon-dominated regions (PDRs) are expected to show a layered structure in molecular abundances and emerging line 

emission, which is sensitive to the physical structure of the region as well as the UV radiation illuminating it. 

Aims. We aim to study this layering in the Orion Bar, a prototypical nearby PDR with a favorable edge-on geometry. 

Methods. We present new maps of 2'x2' fields at 14"-23" resolution toward the Orion Bar in the SO 88-9ci, H2CO 5i5-4i4, "CO 

3-2, C2H 49/2-37/2 and 47/2-35/2, C'*0 2-1 and HCN 3-2 transitions. 

Results. The data reveal a clear chemical stratification pattern. The C2H emission peaks close to the ionization front, followed by 
H2CO and SO, while C"*0, HCN and '^^CO peak deeper into the cloud. A simple PDR model reproduces the observed stratification, 
although the SO emission is predicted to peak much deeper into the cloud than observed while H2CO is predicted to peak closer to the 
ionization front than observed. In addition, the predicted SO abundance is higher than observed while the H2CO abundance is lower 
than observed. 

Conclusions. The discrepancies between the models and observations indicate that more sophisticated models, including production 
of H2CO through grain surface chemistry, are needed to quantitatively match the observations of this region. 

Key words. ISM: molecules - ISM: structure - ISM: individual: Orion Bar - stars: formation 



1. Introduction 

In photon-dominated regions (PDRs), UV radiation between 
a few and 13.6e V drives the thermal and ch emical balance 
of interstellar gas dH oUenbach & TielensI [19991 and references 
therein). PDRs are ubiquitous in the universe: surfaces of molec- 
ular clouds adjacent to Hii regions, planetary nebulae, proto- 
planetary disks and also the nuclei of distant active galaxies. 
As such, PDRs are signposts of radiative feedback processes 
driven by star formation. Heating proceeds through photo- 
electric emission from dust grains and cooling is mostly due to 
fine-structure emissions of [Cii] and [Oi], and CO rotational 
lines (e.g.. [Ste rnberg & Dalgarno 1995; Kaufman etal. 1999; 
lMeiierink&Spaansll2005b iRolliget al.ll2007l) . Shielding of the 
UV radiation by dust and gas creates a layering structure where 
a sequence of different chemical tr ansitions is produced b y the 
gradual attenuation of the UV field dOssenkopf et al.l[2007h . 

The Orion Bar is a prototypical PDR located between 
the Orion molecular cloud and the Hii regio n surroundin; 
the T rapezium stars, at a distance of 414 pc (Menten 



(1 - 4) X 10'*;t'o ( IHogerheiide et al.lll995h . where the Draine 
field ;t'o = 2.7 x 10"^ erg s"' cm"^ 

The dumpiness of the PDR infeiTed by IHogerheiide et al 



(1995|) was confirmed by in terferometric data ( Young Owl et al 
2000riLis & Schilkd |2003|). Cliimp den sities up to 10^ cm " 



iding 
et alj 



2007). Multi-wavelength observations (Tielens et 
1993[ [V an der Werf et al. Il996t i Hogerhei ide et all 119951 



Young Owl et al.u200Ck , Walmslev et al...20 00) indicate a geome- 



try for the Bar where the PDR is wrapped around the H 11 region 
created by the Trapezium stars and changes from a face-on to 
an edge-on geometry where the molecular emissions peak. The 
mean density of the Bar is about 10^ cm"^, the mean molecular 
gas temperature 85 K, and the impinging radiation field is 



were derived by Lis & Schilkel (l2003l) while the density of 
the interclump medium should fall between a few lO'* cm"^ 
(Young Owl et al. 2000) and 2 x 10^ cm^^ (Simon et al. 199j). 
The physical stratification of the PDR is well established 
(iTielens et al.'l993 yVan der Werf et alJl996HSimon et al.ll997l: 
Lis & Schilke 2003|). Vibrationally excited H2 emission is lo- 
cated 15" from the ionization front, where HCO^ 1-0 peaks as 
well, while the CO peak resides at 20" and CS slightly further 
in. 

This paper presents emission maps of various molecular 
species, allowing us to probe and understand the chemical strat- 
ification in the Orion Bar in greater detail. 



2. Observations 

This work is based on data obtained as part of the Spectral 
Legacy Survey (SLS. IPlume et al |2007|) . being conducted at the 
James Clerk Maxwell Telescope[J on Mauna Kea, Hawai'i. The 



' The James Clerk Maxwell Telescope is operated by The Joint 
Astronomy Centre on behalf of the Science and Technology Facilities 
Council of the United Kingdom, the Netherlands Organisation for 
Scientific Research, and the National Research Council of Canada. 
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SLS is performing spectral imaging of 2'x2' fields in the direc- 
tion of the Orion Bar and four other targets. Once completed, 
the spectral range covered by the SLS will be 330-362 GHz. 
This Letter presents the first results from a selection of molec- 
ular lines in the already acquired data: two C2H lines and one 
ti-ansition each of SO, H2CO and '^CO between 330 and 350 
GHz, and two lines in auxiliary data in the 230 GHz window: 
C"^0 2-1 and HCN 3-2 (Table[B. 

The SLS uses the 16-pixel HARP receiver (325-375 GHz) 
and the ACSIS correlator. The angular resolution of the JCMT 
at the observed frequencies is 14-15" (!i!30mpc at the distance 
of the Orion Bar); the harp4_mc jiggle position switch mode 
produces 2'x2' maps sampled every 7'.'5. The spectra are cal- 
ibrated by observations at an off-position 10' to the southeast. 
At an exposure time of 4 minutes per map pixel, the rms noise 
level in the current SLS data is 0.07-0.08 K in velocity bins of 
0.83-0. 89km s ' . Single sideband system temperatures for these 
observations are ~700 K. 

The HCN 3-2 and C'^O 2-1 lines were observed in 2005 
with the A-band (215-275 GHz) receiver and the DAS spec- 
trometer The same 2'x2' field as with HARP was mapped at 
half-beamwidth spacing, using the raster mapping mode. The 
JCMT beam size at 220-266 GHz is 19-23" («;40 mpc). The off- 
position for the A-band observations was at 20' southeast of the 
map center. The HCN and C'^O observations reach a noise level 
of 0.5-0.6 K at a spectral resolution of 0.18 and 0.21 kms"', 
respectively. System temperature values are ~350K (DSB); the 
integration time is 19 seconds per map pixel. The pointing ac- 
curacy was checked every hour and was generally better than 2" 
during all observing runs. 

The data were manually reduced using standard Starlink 
procedures. The SLS observations were handled in blocks of 
1 .6 GHz effective bandwidth and the other data as one set per 
emission line. The time series were converted to data cubes and 
first-order polynomial baselines were subtracted. Line intensi- 
ties quoted here are given in terms of the main-beam tempera- 
ture Jmb, i e-, they have been corrected for the JCMT main-beam 
efficiency of 63% in the 345 GHz window and 69% in the 230 
GHz window. Single-channel signal-to-noise ratios at the peak 
in the northeast sHce are -16 for C"^0 and SO, -40 for HCN, 
H2CO and C2H, and 900 for '^CO. 

3. Results 

Figure [1] presents images of the observed molecular emission, 
overlaid on an 8 fim image of the Orion Bar taken with the 
Spitzer Space Telescope (Megeath et al., in prep.). The mid- 
infrared emission traces warm dust at the surface layer between 
the H II region and the molecular cloud. The molecular emission 
is seen to peak to the southeast of this layer, and the JCMT data 
reveal a clear stratification. The C2H emission peaks close to the 
surface layer, followed by SO and H2CO, whereas CO and HCN 
peak deep into the molecular cloud. 

The stratification is more easily visible in Fig.|2]which shows 
slices through the images in Fig. [1] taken along the arrows in 
Fig- [H perpendicular to the ionization front. To produce these 
intensity profiles, all maps are regridded to 3'.'75 pixels, roughly 
half the pixel size of the original maps. Fitted peak positions and 
peak fluxes for the seven profiles are listed in Table[Tl along with 
the spectroscopic parameters of the observed lines; the position 
uncertainties include Gaussian fitting errors (< 1" » 2 mpc). 

In addition to the stratification, the molecular emission 
shows a two-peak structure. The northeast peak is brighter in 
HCN and H2CO, while the southwest peak is brighter in SO and 
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Fig. 1. Images of molecular emission, observed with the JCMT 
toward the Orion Bar of the species and transitions listed in 
Table[T] Line intensities are integrated over 5-8 km s ' , depend- 
ing on the line width. The C2H transition shown here is the 4ci/2- 
37/2 transition; the other transition shows a similar spatial distri- 
bution. Contour levels are drawn at 10, 20, . . . , 90% of the max- 
imum integrated intensity for every map: 13.3 Kkm s"' for C2H, 
8.7 for H2CO, 7.7 for SO, 19.9 for C'^O, 70.5 for HCN, and 
276.8 for '^CO. Note that the measured maximum J r^bdv val- 
ues across the northeast slices (Table [U are not the same as the 
highest value across the entire map. The grayscale background 
image shows Spitzer Space Telescope 8 fim continuum emission. 
The filled circle in each frame indicates the beam size at the rel- 
evant frequency. Axes are annotated with right ascension and 
declination (J2000). 



the peak brightness is roughly equal for CO and C2H. These 
variations are probably caused by variations in the filling fac- 
tor of dense gas, a s seen, e.g., in the high-resolution data by 
iLis & Schilkd ( 1200 3). The measured line widths of 5-8 km s"' 
will be useful for future detailed modeling, however, this paper 
focuses on integrated intensities. 

We use the obse rved line intensities and the radiative trans- 
fer program RADEX der Tak et al ]|20071 to estimate beam- 
averaged molecular column densities. Although temperatures 
are known to range from 5 K in dense clumps up to 150 K in 
the interclump medium (B atria & WilsonI [20031: iLis & Schilkd 
2003), we assume a kinetic temperature of 85 K and a volume 
density of 10^ cm~^, as a pplicable to the extended molecular gas 
(iHoeerheiide et al.]ll995l) . Since some of the line emission may 
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Fig. 2. Profiles of molecular emission along slices perpendicu- 
lar to the ionization front at position angle 225°(see aiTows in 
Fig. [TJ: (top panel) through the northeast core; (bottom panel) 
through the southwest core. The zero point of the distance 
scale is placed at the C"*0 peak position: (RA=5''35™25'4, 
5=-5°24'37") for the northeast slice and (RA=5'>35™2P.6, 
5=-5°25'19") for the southwest slice. 
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Fig. 3. Top panel: model calculation of relative abundances as a 
function of depth into the cloud. The dashed black curve shows 
the temperature profile. The thin vertical line indicates the zero- 
point of the horizontal scale in the bottom panel. Bottom panel: 
modeled relative intensity profiles in the direction perpendicular 
to the line of sight, convolved by the JCMT beam. In this panel, 
the zeropoint of the depth scale is set to the modeled C'**0 peak, 
at 36 mpc from the edge of the molecular cloud. 



originate from dense clumps, we also list column densities es- 
timated in the high-density (LTE) limit (Table [TJ. Both for the 
10'' cm"-' case as well as the LTE limit, varying the tempera- 
ture by +30 K changes the column densities by <30% (where a 
higher temperature gives a lower column density), so the molec- 
ular column densities are reasonably well constrained. 

Our values for the column densi ties agree with previously 
published values to a facto r of ~2: iHog erheiide et al.l (Il995h 
for H2CO, an d HCN, iJansen etalj ^199 5) for C2H, and 

iLeurini et al.l ( l2006l) for SO. The published column densities 
for HCN and SO are closest to our high-density (LTE) esti- 
mates (Table [TJ; this is consistent with these two species be- 
ing confined to the high-density gas while the H2CO emission 
has contributions from both density phases. The spatial distri- 
bution of the emission in this work is similar to that of pre- 
viously observed single-dish emission from molecules such as 
'^CO, '^CO, HCO+, H? CN and CS (e.g iTauber et alJ[T994l 
IVan der Werf"etan[T996l:[Simon et al.l[l997l) . 

To convert the observed C'^O line intensity to a total 
H7 column dens ity, we assume an isotopic ratio '^O/'^O=500 
dWilson & Roodi ri994) and a CO a bundance of l.lxlO'^ as ap- 
plicable for the Orion Bar PDR dJohnstone et all l2003h . The 



result is A^(H2) = 1.0x10^-' cm ^, which implies the follow- 
ing molecular abundances at the peak locations: x(C2H) = 
2x10-^ -x(H2CO) = 4x10-1"; x{SO) = TxlO""'; x(HCN) > 
5x10"'". Here we have assumed that the emission from all 
m olecules except CO originates in the high-density gas. Based 
on Ijansen et al T (ll995h . who propose the clumps to encompass 
10% of the material along every line of sight, we estimate the 
clumps to provide for roughly half of the total column density 
for those molecules that exist in both phases. Hence, the result- 
ing uncertainty in the derived abundances would not exceed a 
factor 2. Since C2H and CO are hardly coexistent, the C2H abun- 
dance is probably underestimated by a factor of 3-10. 

4. Discussion 

To interpret our observations we compare them to the pre- 
dictions from a PDR model. Previous models of the Orion 
Bar do not present abundances for all six species in our data, 
so we use the global conditions of the Bar {x-3xlO'^xo^ 
mh^IO^ cm"^, see Sect. [TJ to construct a simple model us- 
ing the KOSMA-T PDR code (iRoUig et al.1 l2006h . We refer 
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Table 1. Observed molecular lines. 



Molecule 


Transition 


Frequency 




Peak position (mpc)''*' 


/r,„bdv<'=' A'„,b (cm-2)<i5) 






(MHz) 


(K) 


NE slice SW slice model 


(Kkms-') n=10^cm-^^ LTE 





Nj 


= '^7/2 


— > 35/2 


349400.5'"> 


42 


-14 ± 


1.8 


-23 ± 1.0 


-16 


8.6 ± 1.6 


4.4x10"* 


2.3x10'-* 


C2H 


Nj 


= 49/2 


— > 37/2 


349337.8*'" 


42 


-14 ± 


1.5 


-23 ± 1.6 


-16 


11.3 ± 2.0 


4.9x10"' 


2.4x10'" 


0-H2CO 


J KpKo 


= 5,5 


^4i4 


351768.6 


62 


-3 ± 


1.3 


-7 ± 1.9 


-12 


9.3 ± 1.6 


7.7x10'^ 


3.8x10'^ 


SO 


Nj 


= 89 




346528.5 


79 


-2 ± 


1.4 


-5 ± 1.0 


8 


5.2 ± 0.9 


1.4x10'^ 


7.2x10'^ 


c"*o 


J 


= 2 


^ 1 


219560.4 


16 


0± 


1.2 


+ 1.4 





20.3 ± 3.4 


2.3x10'* 


2.9x10'* 


HCN 


J 


= 3 


2 


265886.2*'^' 


26 


3 ± 


2.5 


-1 ± 1.8 


8 


74.0 ± 14 


>6.5xlO'2<f 


>4.9xl0'^<f 


'^CO 


J 


= 3 


^ 2 


330588.0 


32 


4± 


1.9 


-1 +3.0 


2 


186 ±32 


>2.1xl0'*<" 


>2.5xlO'*<''' 



(a): Blend of F=4-3 and F=3-2 hyperfine components (separation 1.4 MHz), (b): Blend of F=5-4 and F=4-3 hyperfine components (separation 
1.3 MHz), (c): Blend of F=3-3, F=3-2 and F=2-l hyperfine components (total separation 1.5MHz). (d): Fitted position of the emission peak 
along the slice through the northeast and southwest cores (see Fig.[T](, and the modeled slice (Fig.[3] bottom panel), all relative to the C'^0 peak, 
(e): Fitted value of integrated intensity at peak position along the northeast slice (see Fig.[T]l; uncertainties include fitting uncertainty and 15% 
absolute calibration uncertainty, (f): Column density uncertain due to high optical depth; values listed here are lower limits, (g): Beam-averaged 
column densities for n=lO^ cnT^ and in LTE, both at 7=85 K. 



to iRolli^ et"an (l2007h for a comparative study of other PDR 
codes. To simulate the plane-parallel geometry of the Orion Bar 
we represent it by one side of a very large sphere. Starting 
from atomic abundances of ^(0)^3.0x10 '"', jc(C)=1.4xlO~"' 
4S)=2.8xlO"^ and jc(N)=1.0xlO (see also lCubick et alJ2008l 
Table 3), the model computes the equilibrium chemistry and 
temperature of the clou d usin g reaction rates from the UMIST 
database (Woodall et al.l l2007h . Based on the geometry in Fig. 1 
of iJansen et al.i (.1995i) we simulate the PDR as a layer of thick- 
ness 49 mpc where the Bar represents a part of the layer that 
is parallel to the line of sight, so that its apparent depth is in- 
creased by a factor 10 compared to its thickness. We note that 
the critical densities are ~10^ cm"^ for all transitions in Table [1] 
except thos e of th e CO isotopes, while the highest H2 density 
reported bv lLis & Schilke (2003) is also 10^ cm"^. Hence, to ap- 
proximate the spatial distribution of the line emission, we as- 
sume optically thin LTE and weight the abundance profiles with 
a Boltzmann cutoff: /rei - (x/xq) exp(l - Tgas/^up) if T'gas < E^p 
and /rei - x/x() otherwise. Here /jei is a normalized intensity, x is 
abundance, Tgas is the modeled local gas temperature (see Fig.|3] 
top panel), E^p is given in Table [1] and the normalization factor 
xq is chosen from line to line. The distribution is convolved with 
the JCMT beam at the appropriate frequencies to obtain the pre- 
dicted relative intensity profiles (Fig. [3] bottom panel). 

Table[T]compares modeled positions of the peak emission to 
fitted peak positions in the observed profiles for each species. 
Overall our simple PDR model matches the observed layering 
of C2H-H2CO-CO. For H2CO the highest abundance occurs in 
the cold shielded material, but the warm outer layers dominate 
the emission due to the higher £up for the observed transition. 

There are also significant discrepancies between the model 
and the observations, (i) The model predicts that SO peaks deep 
in the cold material (at 8 mpc), following the distribution of 
HCN, not at the location where H2CO peaks (in front of CO) 
as is observed, (ii) Conversely, the H2CO peak in both slices is 
observed closer to the peaks of the CO isotopes than to the C2H 
peak, while the model predicts H2CO to peak significantly closer 
to C2H. (iii) For HCN the model peaks 5-9 mpc deeper than is 
observed and (iv) the observations show an HCN emission wing 
toward the PDR surface which is weaker in the model, (v) We 
find that the absolute abundance of H2CO and the lower limit to 
the absolute abundance of HCN at the model peak are consis- 
tent with the observations, while the model abundance of C2H 
is ~2 times smaller than observed, (vi) In addition, our modeled 



SO abundance is orders of magnitude larger than the observed 
value. 

The underestimate of the C2H abundance predicted by 
our pure gas-phase model might be explained by C2H be- 
ing forme d via photo-destr uction of PAH molecules, as pro- 
posed by iPetv et al] (l2005h . This explains a detectable C2H 
abundance in the outer la yers where the PAH emission peaks 
(IVan der Werf et al.l[T996h . 

Our modeled SO abundance is a factor 1.4x10"^ larger than 
observed. This is explained by the too high gas phase sulfur 
abundance of 2.8x10"^ used in the KOSMA-t model, whi ch rep - 
resents diffuse clouds and no depletion ( Federman et al.,ll993h . 
Recent estimates of cosmic sulfur abu ndance provide values as 
low as 1.4xl0~^ dAsplund et al.ll2005l) and sulfur depletion fac- 
tors in de nse clouds may rang e from 4 dOoicoechea et al. 200^ 
to 1000 ( Tieftrunk et al.lll994i) . Depletion factors of several 10^ 
would be required to match the Orion Bar observations with the 
current model. The low depletion seems not to apply to the entire 
Orion region, but to be specific for the Horsehead, which has a 
particularly low density. 

The abundance of H2CO derived in the LTE limit is slightly 
higher than the modeled abundance; the discrepancy increases 
sharply as soon as the H2CO gas is not completely in LTE. This 
abundance discrepancy might be explained by the evaporation 
of molecu les from warm ( ~30K ) dust grains. Based on reaction 
rates from I Woodall et aU (|2007|) . the photo-dissociation lifetime 
for H2CO at Ay - 5 mag is only ~10yr, which introduces the 
need for a continu ous supply of H2CO, e.g., from grain sur- 
faces (Herbst 2000). This is consistent with the suggestion by 
Leurini et al. ( 2006) that H2CO may originate primarily from the 
warm interclump medium. 

Finally, the observed broad spatial profiles for HCN and 
'^CO (Fig. |2]l may be explained by optical depth effects, since 
the optical depth inferred by RADEX is 5-10 for these molecules 
and <1 for the others. The different spatial distributions of the 
two CO isotopologs result mainly from the different beam sizes. 

Future observations will test our hypothesis that grain sur- 
face chemistry plays a key role in PDRs. One such test will 
be JCMT-SL S obse rvations of CH3OH in the Orion Bar: while 
iLeurini et al.l (l2006l) suggest that CH3OH traces the clumps and 
H2CO traces the interclump medium, both molec ules are thought 
to be the products of CO hydrogenation on ice (Watanabe et alj 
2004), so their spatial coincidence would mark a critical vali- 
dation step. Other analysis will explore sulfur chemistry using 
sulfur-bearing species covered by the SLS, such as CS and SO2. 



M. H. D. van der Wiel et al.: Chemical stratification in the Orion Bar: JCMT-SLS data 



5 



Once completed, the SLS data set will enable studies of mul- 
tiple transitions of various molecules, so that the excitation of 
molecules (as well as their velocity structure) can also be com- 
pared to models. In the longer term, (sub)millimeter-wave imag- 
ing spectroscopy of other Galactic PDRs is needed to establish 
how common grain chemistry is and how it depends on envi- 
rormient. Important tests wiU be spectral surveys with the HIFI 
instrument on board ESA's Herschel Space Observatory, which 
will target hundreds of molecular lines in many Galactic PDRs 
and allow a detailed picture of their chemistry to emerge. 

Acknowledgements. The authors acknowledge the JCMT staff for their support, 
Markus RoUig for help with the PDR model, and Peter Schilke, Maryvonne 
Gerin, and an anonymous referee for their careful reading of the manuscript. 
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